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Brane- World models suggest modification of Newton's law of gravity on the 3-brane at submil- 
limeter scales. The brane-world induced corrections are in higher powers of inverse distance and 
appear as additional terms with the Newtonian potential. The average inter-particle distance in 
white dwarf and neutron stars are 10~ 10 cms and 10 -13 cms respectively, and therefore, the effect 
of submillimeter corrections needs to be investigated. We show, by carrying out simple manybody 
calculations, that the mass and mass-radius relationship of the white dwarf and neutron stars are 
not effected by submillimeter corrections. However, our analysis shows that the correction terms in 
' the effective theory give rise to force akin to surface tension in normal liquids. 
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I. INTRODUCTION 



String theory is considered as a strong candidate for the unified theory of the standard model and gravity. However, 
string theory is a higher dimensional theory, and therefore, it is necessary to find a mechanism to get rid of the extra 
C*~) • dimensions with some imprint left on the four dimensional physical world. But such effects must be consistent with the 
existing experimental and observational data. The brane-world scenario proposed by Randall and Sundrum0, 0] for 
' describing a low dimensional physical theory starting from a higher dimensional theory has attracted lot of attention 

recently The model has triggered a large amount of work in cosmology and astrophysics Q. 
, The brane-world scenario leads to modification of gravity at submillimeter distance scales. The brane corrections 
are in higher powers of inverse distance and appear as additional terms in the Newtonian potential. In the Randall- 
, Sundrum scenario, the modification suggested for the Newtonian potential is given by, 



O" 1 r 3r 2 



$( r ) = (1 + — j) (1) 



where I is the curvature scale of 5-dimensional anti-deSitter space time, and G is the usual Newton's constant of 
gravity, M is the mass, r is the distance in 3-space. In the setting above, I has the dimension of length in 3-space, in 
other words, it is measured in centimeters in standard CGS units. Henceforth, we will take l 2 s = |i 2 

Several groups have estimated bounds on the distance scales at which such corrections may be allowed without 
violating the existing experimental and observational data Q . As mentioned above, there have also been some works 
dealing with the astrophysical consequences of such modifications of gravity 0, S 0- The average inter-particle 
distance in white dwarf and neutron stars are I0 -10 cms and 10 -13 cms respectively, and therefore, the effect of 
submillimeter corrections needs to be investigated. Wiseman has carried out a detailed numerical study of the 
neutron stars on the brane and claimed that the normal parameters of the neutron stars such as mass and radius 
are not effected by brane world corrections. In this paper, we carry out simple manybody analytical calculations 
using the brane modified gravitational potential. Our approach is based on manybody theory in semi-relativistic 
quantum mechanics [jj, Il0|. At first, we briefly summerise the nature of the potential and the force that emerges in 
the manybody context due to the barncworld correction to gravity (the detail calculation is presented in appcndix-I). 
We find that, for extended massive bodies, the nature of the force and the potential due to brane correction terms 
is very different from those due to Newtonian gravity. For example, inside a hollow spherical shell the potential due 
to Newtonian gravity is constant, independent of location of the test particle. Therefore, the force is zero inside the 
shell. We show that this is not true in the presence of the brane correction term. This result may have some relevance 
in context of experimental observation of deviations from Newton's law of gravity. We also show that the brane-world 
correction to gravity gives rise to force akin to surface tension in normal liquids. Following that we consider the effects 
due to the brane-world corrections on the mass and the mass-radius relationship of white dwarf and neutron stars. 
We find that the brane world corrections do not influence these two parameters in any significant way. 
In the appendix-I, we provide detailed calculations of the potential, potential energy and the force due to the brane 
correction term. We also point out that many new physical features that emerge from corrections to Newtonian 
gravity are generic due to the deviation from the inverse square law of force and does not depend on the particular 
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form of the correction term. 

In appcndix-II, we analyse the nature of the extrema of various energy functionals that appear in the text. 



II. MANY-BODY EFFECTS OF BRANE WORLD CORRECTION TO GRAVITY 

We consider in this section the many-body effects of brane world correction to gravity in a classical theory. Essen- 
tially, we carry out analysis of the potential and force due to brane world correction for spherically symmetric body of 
constant density and a fixed radius. Details of calculations are given in appcndix-I. Here we provide a brief summary 
of these results. As already mentioned in the introductory section, for extended massive bodies, the nature of the 
force and the potential due to brane correction terms is very different from those due to Newtonian gravity. These 
claims are demonstrated latter in this section. 

For a spherically symmetric ball of constant density p and radius a, the potential due to brane corrections at a point 
f ( r > a , r being the distance from centre of the ball) outside the ball is given by, 

-r i \ _ „ , 2 r , r + a 2a 

$bW = -2TxGlip\ln — 

r - a r 

The corresponding force is given by, 

F = 2irGl 2 s p [— — + H£] = -AnaGll p [-^—^ - \] 

s r + a r-a r 2 3 r 2 - a 2 r 2 ' 

When the point r is close to surface of the ball, both potential and the force grows and becomes very large. Therefore, 
one needs to introduce a parameter e which characterises the smallest inter-particle distance. At a distance e from 
the surface of the ball, the force is given by, 

F = 2nGl 2 s p [-- + -] 
e a 1 

This is clearly a skin effect. However, this force is not as strong as one may think. At a distance of 1 mm from the 
surface of a neutron star, this force is weaker than the force of Newtonian gravity. The is because this force depends 
on the density and not on the mass. However, this force can be large close to the surface of a spherical grain of very 
small radius and may even overtake the force due to Newtonian gravity. 

For a spherically symmetric shell of outer radius R and inner radius b, the potential at an interior point f (r < b ) 
is given by, 

d>' s (r) = -2nGl 2 p In^r^ 

Force due to this potential, at distance r from the centre is given by, 

9r 9r 
F{r) = 2,Gllp [— — - -—-] 



From the expression above it is obvious that F(r) > 0. This means that the force is directed in the direction in 
which r increases. In other words, it is directed towards the inner boundary surface of the shell. At the centre of 
the shell, r = 0, and therefore, the force is zero. But any small perturbation will move the test particle away from 
the centre towards inner boundary surface. This is a new element that emerges from the brane-world correction to 
Newton's law of gravity. We know that the force due to Newtonian gravity is zero inside a spherically symmetric shell. 
If one can find the above mentioned effect inside a shell, it will provide a very clean test of brane world corrections to 
Newtonian gravity. May be such experiments can be carried out on board artificial satellites. 

In this case, as in the case of spherical ball, there is skin effect. The force of attraction close to the inner boundary 
of the shell at distance, say, e from the inner boundary surface is given by, 

F = 2nGl 2 p [\ jjgL-] 

Therefore, any interior layer of a spherical ball experiences two types of forces: one is due to a sphere just bellow 
it and other is due to the spherical shell outside. These two types of forces, individually, give rise to a large forces 
on this layer due to the skin effect. Remarkably, these two forces act in opposite directions and the forces due to 
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skin effect are equal in magnitude and exactly cancel each other. On the layer at r, let F\ be the force due the inner 
sphere and F 2 be the force due to the outer shell. 

F 1 = 2-kG1 2 s p [-- + -] 
e r 

1 9r 

Therefore, the net force at any point on the layer at a distance r from the centre of the ball is given, 

F ~ Fi + F 2 = 2nGl 2 s p [- - ~^-^\ 
r R z — r l 

This result is valid for every layer except the one at the outer most boundary. For the outer most layer, there is an 
inner sphere but no outer shell. There is no cancellation of the skin effect on this layer. Therefore, it experiences a 
large force. This will give rise to effect similar to surface tension in normal liquids. 

The formula above also shows that the net force due to the brane-world induced corrections at any point on layers 
at distance r where r < is repulsive, and on layers at a distance r , r > is attractive. Within a submillimeter 

central core of a neutron star there are approximately 10 36 neutrons arranged in approximately 10 12 layers. At 
any point on these layers the net repulsive force due to the brane-world induced corrections overtake the attractive 
force due to Newtonian gravity. Also note that any change in the radius of the star due to mass acreation or mass 
loss leads to the rearrangement of the field of force in the interior of the star. Therefore, the interior and the exterior 
of the relativistic stars are intimately (although very weakly) connected to each other in the presence of brane-world 
corrections. 

The potential energy of a spherically symmetric body (of constant density) of mass M and radius R, due to the 
brane-world induced corrections is given by, 

3GM 2 l 2 R 3GN 2 m 2 l 2 sl R 

7 = R^ ln l 

Here N is the number of particles in the spherical ball of radius R, m is the mass of the elementary constituent 
particles, In a classical theory, the parameter e should be considered as the minimum inter-particle distance. But in 
a quantum theory, as we will show in the latter section, it is possible to fix its' value using Heisenberg's uncertainity 
relation and Pauli's exclusion principle. 



III. MASS-RADIUS RELATIONSHIP OF THE WHITE DWARF AND NEUTRON STARS ON THE 

BRANE 

In this section, we show that the brane world corrections to Newton's law of gravity does not influence the 
mass-radius relationship of white dwarf and neutron statrs in any significant way. Let us consider a system of large 
number of fermions contained in some volume V. Let us also assume that the system is almost a compact sphere 
with a large constant density such that fermions are relativistic. We will consider both the white dwarf and neutron 
stars. We do not discuss here the details of the nuclei present in the white dwarf stars. We simply consider that, in 
the white dwarf stars, there are two types of particles: the electrons and the nucleons. For the nucleons, we have in 
mind the protons and the neutrons. The over all charge neutrality of the star implies that number of electrons and 
the protons should be equal. This suggests that for N electrons in the star there should be 27V nucleons ( N protons 
and N neutrons). However, we will consider that the number of electrons is equal to the number of nucleons which 
are essentially protons. Such a choice changes the mass and radius of the white dwarf star by a small numerical 
factor but the analysis is easily carried over to the case of neutron stars. The zero-point kinetic energy in the white 
dwarfs is solely due to the zero-point motion of the electrons but the gravitational potential energy is soley due to 
the gravitational interaction of the nucleons among themselves. For the neutron stars the kinetic as well as potential 
energy is generated by the neutrons alone. We assume that the mass of a nucleons is equal to the mass of a neutron. 
We also consider a situation when the kinetic energy due to temperature is much small compared to the zero point 
kinetic energy and, therefore, can be neglected. 

To make our derivation more transparent, at first, we derive the mass-radius relationship without the brane cor- 
rection terms. Let us consider first the white dwarfs. We have a spherically symmetric region of radius R containing 
N electrons and N nucleons. We assume the electrons to be relativistic; their total kinetic energy is, 

K N = N c(p 2 + m 2 e c 2 ) 1/2 - Nm e c 2 (2) 
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Substituting the value of zero-point momentum, p = hN R / in the expression of the kinetic energy, we obtain the total 
energy as 

™™ »r ^ 2 N 2/3 2 2 , 1/2 A7 2 3 Gm^iV 2 

== AT C (—^2— + m e C ) ' - Nm eC - g X 1 

The quantum mechanical variational ground state is obtained by demanding, dE d ^ — O.This gives, 

N^h 3 GmlN 2 /* 
= — x — (3) 

(h 2 N 2 / 3 + m 2 e c 2 R 2 y /2 5 ch K ' 

Both the LHS and RHS of the equation above are positive. The LHS, as is clear from equation above, is strictly less 
than one, and therefore, the RHS is strictly less than one. In particular 

3 GmlN 2 ^ 
5 ch 

This provides an upper bound on the number of particles and therefore, an upper bound on the mass of the star. As 
before, the upper bound on the number of particles, N max , is given by, 

„ r9 /o 5 ch 
N ' = - x 



3 Gm 2 n 

Therefore, as before the upper bound on the mass of the star is M max — m n N max - Let us call the radius defined as, 

R — 



m e c 

Compton radius of a compact relativistic white dwarf. Let us also introduce the notation, 

a= {N/N max ) 2 ^ 

Using this, we can write Eq.(19) as 



(i+ %Y' 2 



This can also be expressed as, 



R = R \j\-l (4) 



This is the mass-radius relationship of white dwarfs in the absence of brane correction term. (It is shown in 
appendix-II, that this solution corresponds to the minimum of the total energy, as a function of R). Note that a, by 
definition, is strictly less than one. We can get another bound on a as follows. The white dwarf star that we are 
considering is relativistic, and therefore, p = — > m e c. This implies that R < h ^ c = R a . Taking R < R , 
in the equation above, we find that a > An order of magnitude calculation gives the mass of the white dwarf 
between one to two solar masses and the radius a few thousand kilometers. Neutron star parameters are obtained 
simply by replacing the electron mass by neutron mass. As an order of magnitude estimate, the mass of the neutron 
stars turns out to be similar to the white dwarf mass but the radius is just a few kilometers. 

Let us now consider the effect of the brane-world induced correction term on the parameters of white dwarf and 
the neutron stars. 

The potential energy of spherically symmetric ball due to Randall-Sundrum brane corrections is given by, 



3Gls 2 m 2 n N 2 . R 
R? 

Therefore, the total energy including brane correction is given by, 



Urs = x In- 

R A e 



E(R) =K N + U N (R) + U RS = N c (p 2 + m \c 2 fl 2 Nm e c 2 | x _ ZGls 2 m 2 n N 2 x ^R_ 



5 



Substituting p — hN R /3 , we obtain 

»r ,K 2 N 2/3 2 2Nl/2 A7 2 3 Gm 2 n N 2 3Gls 2 mlN 2 , i? 
E(R) =N c n9 + mlc 2 ) 1 ' 2 - Nm e c 2 - - x 2 _2 x Zn- 

The quantum mechanical variational ground state is obtained as before, 

NV 3 h 3 GmlN 2 / 3 9Gls 2 m 2 n N 2 / 3 R 

x IS + 5di >< (5) 



(h 2 N 2 / 3 + m 2 e c 2 R 2 ) 1 / 2 5 c?i i? 2 cfr e 

Note that e was introduced as shortest inter layer distance in the spherical ball. From physics point of view it can 
be identified with the shortest inter-particle distance. For fermions, as per Pauli's exclusion principle, it can not be 
shorter than R/N 1 / 3 . We take e = R/N 1 / 3 ( In f = InN 1 / 3 ). 

Both the LHS and RHS of the equation above are positive. The LHS, as is clear from equation above, is strictly 
less than one, and therefore, the RHS is strictly less than one. In particular 

3 x GmlN 2 / 3 { 
5 ch 

This provides an upper bound on the number of particles N, and therefore, an upper bound on the mass of the star, 
M = m n N. We will proceed in the same way as we did in the absence of the brane correction terms. However, in 
this case we obtain a cubic equation for the square of the radius. Before we solve the equation, let us introduce the 
following notation, 

n 15ls 2 , R Ihls 2 , , rl/ , 

f3 = — =- x In— = — =- x InN 1 ' 3 (6) 



We introduce a new variable, 



R 2 e R 2 



Y=A 2 



The equation Eq. J5j), in terms of the new variable can be written as, 

1 _ a(3 

(i + y)V2 - a + ^r 

After some algebraic manipulations, we obtain 

Y 3 + (2/3 + 1 - \)Y 2 + (2/3 + (3 2 )Y + (3 2 = (7) 

Before we look for solution of this equation, let us have, an order of magnitude, estimate of the values of parameters 
a and [3. They are both pure numbers. As mentioned earlier, both for white dwarf and neutron stars, the range of 
a is given by < a < 1. The parameter /3, depends on the Compton radius of the star, and therefore, its' value is 

different for white dwarf and the neutron stars. Taking N « 10 57 , we find that for white dwarfs, (3 w 10 -16 , and for 
neutron stars, (3 w 10~ 10 . The real root of Eq.Q leads to the following, 



Ro V a 2 



H ' ' 1 + 0(0), (8) 

o 



To make the dependence on the brane world parameter l s explicit, we can also write the equation above as, 



These equations are essentially same as Eq.(4) because corrections are insignificant. Note that Eq.(7), is a cubic 
equation, and therefore, there are two more solutions of the equation. Up to order /3, these two solutions are 
coincident, and are given by 



Ho H 
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It is shown, in appendix-II, that these solutions correspond to the maxima of the total energy as a function of R. 
They clearly do not lead to any stable configuration of the compact stars. As shown in appendix-II, the solution 
given by Eq.(8), corresponds to the minimum of the total energy as a function of R. Therefore, in the effective 
four dimensional theory, the standard parameters of the white dwarf and neutron stars do not depend on the brane 
world induced correction term. However, we note that the brane world corrections give rise to surface tension which is 
absent in Newtonian gravity as well as general theory of relativity. This surface tension is the new element in gravity. 
However, we still do not know whether it has any observable signature. 



IV. CONCLUSIONS 



In this paper we have examined the influence of brane world corrections on the parameters of white dwarfs and 
neutron stars using manybody theory. The brane corrections to Newtonian potential give rise to interesting new 
effects, however, these corrections do not influence the stability of the stars in a significant way. For instance, the 
mass, radius and the mass-radius relation remains unchanged. We have demonstrated these features using a simple 
analytical treatment. Our results are in qualitative agreement with the numerical work by Wiseman 0. 

We also point out many new physical features that emerge from corrections to Newtonian gravity. These physical 
features are generic due to the deviation from the inverse square law of force and does not depend on the particular 
form of the correction term: (a.) In the many body context, force as well as potential due to the brane correction term 
depends on density rather than the mass. Therefore, we expect that the relative strength of the force compared to 
Newtonian force will be more close to the surface of spherical grains of smaller radius, (b.) Newton-Birkhoff theorem 
is violated by the brane world induced correction terms. Therefore, inside a hollow spherical shell there will be force. 
This may turn out to be useful for experimental verification of Randall-Sundrum or any other similar corrections to 
gravity, (c) Violation of Newton-Birkhoff theorem due to the presence of brane-world correction terms, gives rise to 
force similar to surface tension in normal liquids. Neither Newtonian gravity nor the general theory of relativity give 
rise to force akin to surface tension. The correction term shares many similarities with molecular force in normal 
Liquids. At present, we are not very clear which of the new features discussed here are relevant to observations. In 
our opinion, this is an important issue which requires further investigation. 
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VI. APPENDIX-I 



We want to calculate the gravitational potential energy of a spherically symmetric ball of radius R and constant 
density p. Here the two-body gravitational potential energy consists of two parts. First part is the standard Newtonian 
part given by, 

U N = (9) 
\r — r | 

where r and r' are the locations of the point masses m\ and mi respectively. The Newtonian potential at a point r 
due to a point mass m at point r' is given by 

*N(?) = -T^r (10) 
\r — r | 

The second part of the potential energy comes from the brane world correction to gravity which in the Randall- 
Sundrum scenario is given by, 

U B = - Gl >l m 2 (11) 
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The brane world correction to the potential at point f clue to a point mass m at the point r' is 

Gl 2 m , . 

= -—^r- (12) 

\r — r\ 6 

For a mass distribution with density p(r'), the formulas above take the following forms, 

p(r')d 3 r' 



J \r — r' 



®N(r) = -G (13) 



<S> B (r) = -Gl 2 j*WL (14) 
' 1 |F-r'| 3 



Un = _ g , , p(rlp(P)dWr> = , $w(f)p(f) ^ (15) 



Us = / / = / »s W/KO^ (16) 

A. Newtonian Potential and Potential Energy for a spherically symmetric ball of constant density 

Let us consider a spherically symmetric ball of radius R and of constant density p. In spherical co-ordinate system 
the formula for $_/v(r) gives, 

_. . , _ f f f r' 2 dr'sin9d9d<j) 

<Mr) = -GpJ J J (r2+r , 2 _ 2rrW)1/2 (17) 

In the expression above, the integral over (f> ranges from to 2n and for 6, it ranges from to w. The integral over <p 
gives a factor of 2-k. The integral over 9 can be carried out by making the following change of variable 

£ = r 2 + r' 2 - 2rr'cos9; d£ = 2rr' sin9d9 ; 

The integration over £ ranges from £ = (r — r') 2 to £ = (r + r') 2 After carrying out integrations over and 9, 
we obtain 

$jv(r) = -2irGp J[{r + r')-\r-r'\] r -dr' (18) 

For r > a (> r'), and < r' < a 

$jv(r) = -2-xGp \ \{r + r') - (r - r'))—dr' = -4w— r u dr' = — = (19) 

Jo r r J 3r r 

For r < 6 (< r'), and b<r' < R 

$jv(r) = -2nGp [ {(r + r') - (r 1 - r)]-dr' = -4irGp [ r'dr' = -2nGp(R 2 - b 2 ) (20) 
Jo r A 

which is constant independent of r . Since force is given by the gradient of the potential with respect to r, F = 
—grad r [<&jv (r)], there is no gravitational force within the spherical shell. Once we know the potential , it is easy to 
calculate the potential energy. It is given by, 

p{r)p{r l )d?rd?r 
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Let us consider the case when r > r' . Taking a = r — e, we obtain 

= ^ — - 

The potential energy is given by, 



U N == [ <S> N (r)p(r)d 3 r 

J(\r\<R) 

Substituting the expression for <&Ar(r), we obtain, 

U N = - — — x Airpr 2 dr 

This integral above is convergent, and therefore, we can take the limit e — > before evaluating the integral. After 
integrating the resulting expression, we obtain 

16tt 2 G P 2 R 5 3 GM 2 , , 

Un = li = ~5 * ~R~ (21) 

B. Potential and Potential Energy for a spherically symmetric ball of constant density due to 
Randall-Sundrum Brane world correction to Gravity 

The potential and the potential energy due to the brane world correction term is calculated exactly in the same 
way as in the case of Newtonian gravity.However, in this case there is a new element. The outer shell exerts force on 
the inner shells, and therefore, the effect of outer shell on the inner shell can not be neglected. 

1. The Potential at a distance from a spherically symmetric ball 

The potential $s(r) in spherical co-ordinate system is given by, 

. . „~ f f f r' 2 dr'sin6d6d<j) 

As in the previous section, the integral over <f> ranges from to 2tt and the integral over 9 ranges from to 7r. Integral 
over (j> gives factor of 2tt. Integral over 8 can be carried out by making the following change of variable 

£ = r 2 + r' 2 - 2rr'cos8; d£ = 2rr'sin8d8 ; 

The integration over £ ranges from £ = (r — r') 2 to £ = (r + r') 2 After carrying out integrations over cj> and 
6, we obtain 



$ B (r) = - 2nGl 2 p ( [— ^- l —]-dr' 

8 J l \r~r'\ r + r'r 

For r > a (> r'), and < r' < a 



(23) 



4> B (r)= -2nGl 2 sP l [-J-^ - -L^f-dr' (24) 







After some algebraic manipulations, this integral can be evaluated. The result is 

$ B (r)= -2TiGl 2 p\ln 25 

r — a r 

Large distance limit of the potential can be evaluated in two different ways. From the formula, 

ra p{r')d 3 P 



t> B (r) = -Gl 2 f 
Jo 
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we find that for r >> a, 

* Gil f a ^Gl 2 pa 3 Gl 2 s M 

The other way is to consider the formula, 

$ B r = -2irGlip[ln - — 

r - a r 

which has been derived for r > a, and therefore, can also be written as, 

$u(r) = - 2ttG I 2 p [ ln 1 + a ^ - — 1 = - 2ttG J? p [ Jn(l + o/r) - Jn(l - a/r) - — ] 

1 — a/r r r 

Replacing the logarithms above by their series expansion, after some algebraic manipulations we obtain, 

4irGl 2 pa 3 4irGl 2 ,pa 5 AnGl 2 pa 7 Gl 2 M 
$ R M = — — — = 



2. Force at a distance from a spherically symmetric ball 



Force at a distance from a spherically symmetric ball can be calculated by evaluating the gradient of the potential, 
F = -Grad^ B {r) = -^g^. The potential is given by, 

$ B (r = -2nGl 2 s p[ln 

r — a r 

and therefore, the force is given by 

F = 2ttG1 2 p I— — + H£] = -AnaGl 2 p [-^-^ - \] (26) 

s r + a r-a r 2 1 sr r 2 -a 2 r 2 ' 

Let us calculate the force at a very small distance from the surface. We take r — a = e, e very small, and substitute 
it in the equation above. After some algebraic manipulation we obtain, 

F = 2nGl 2 p [~ + | (27) 

This is clearly a skin effect. However, this force is not as strong as one may think. At a distance of 1 mm from 
the surface of a neutron star, this force is weaker than the force of Newtonian gravity. The reason is that this force 
depends on the density and not the mass. However, this force can be large close to the surface of a small spherical 
grain and may even overtake the force due to Newtonian gravity. 



3. Potential inside a Spherical Shell 



Let us consider a spherical shell of outer radius R and inner radius b. At a point r inside the shell, the potential is 

$' B (r) = - 2nGl 2 p At^t, - -^ 1 ]-dr' 
J b \r — r\ r + r r 

Since r < b < r' < i?, the formula above becomes, 

$' B (r) = -2nGl 2 sP / [- —^}-dr' - -2nGl 2 p / ^ = -2nGl 2 p In- (28) 

Jb r — r r + r r J b r 2 — r A b 2 - r 2 
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4- Force inside the shell 

Force inside the shell can be calculated by taking the gradient of the potential, F = —Grad [$' B (r)] = — — 
For the potential, 

*' B (r) = -2nGl*p ln^£ 

we obtain, 

F (r) = 2nG? s p [-^ ^-,] (29) 

From the expression above it is obvious that F(r) > 0. This means that it is directed in the direction in which r 
increases. In other words, it is directed towards the inner boundary surface of the shell. At the centre of the shell, 
r = 0, and therefore, the force is zero. But any small perturbation will move the test particle away towards the inner 
boundary surface of the shell. This is a new element that emerges from the brane world correction to Newton's law 
of gravity. We know that the force due to Newtonian gravity is zero inside a spherically symmetric shell. If one can 
find the above mentioned effect inside a shell, it will provide a very clean test of brane world corrections to Newtonian 
gravity. May be such experiments can be carried out on board artificial satellites. 

In this case, as in the case of spherical ball, there is skin effect. The force of attraction close to the inner boundary 
of the shell becomes large. This can be seen as follows. Let b = r + e . This simply means that the layer at r is just 
bellow the inner boundary of the shell (e distance away from the inner surface) . 

In this case the force at r is given by, 

1 2r 

F = 2nGf s p[-- w —^] (30) 



5. Cancellation of the Skin Effect on the Inner Layers of a spherical solid ball 

Any interior layer of the spherical ball experiences two types of forces: one is due to a sphere just bellow it and 
other is due to the spherical shell outside. These two types of forces give rise to a large force on this layer due to 
the skin effect. Remarkably, these two forces act in opposite directions and the forces due to skin effect are equal in 
magnitude and exactly cancel each other. It is easy to verify this. On the layer at r, let F\ be the force due the inner 
sphere and F 2 be the force due to the outer shell. 

F 1 = 2ttG1 2 sP [-- + -] (31) 



F 2 = 2*Gllp [- W — I ] (32) 

Therefore, the net force at any point on a layer at r, is given by 

2 2r 

F = F 1+ F 2 = 2^Gl\ P [- - w —z] (33) 

This result is valid for every layer except the one at the outer most boundary. For the outer most layer, there is an 
inner sphere but no outer shell. There is no cancellation of the skin effect on this layer. Therefore, it experiences a 
large force. This will give rise to effect similar to surface tension in normal liquids. 



6. The Potential Energy 

We found in the previous sections that, for a sphere ball, the potential on a given layer at a radial distance r from 
the centre consists of two parts. First part is due to a interior spherical ball. It is given by, 

2 r , r + a 2a 



$u(r)= -2wGljp[ln 



r — a r 
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The second part is due to the outer spherical shell. It is given by, 

R 2 



$' B (r) = -2nGlipln- 



b 2 — r 2 

The potential energy also consists of two parts. The first part of the potential energy is obtained by taking the radius 
of inner sphere a = r — e in the formula for the potential, multiplying the potential by Airr 2 drp, the mass of the thin 
layer and then integrating the resulting expression from r = e to r = R, where e is the smallest inter particle distance, 
R is the total radius of the spherical ball. This is given by, 

U 1= -2nGl 2 sP l R \ln^—^- 2 ^~ e \ ^pr 2 dr (34) 

J e ^ T 

The second part of the potential energy is obtained by taking the inner radius of the spherical shell, b — r + e in the 
formula for the potential due to the shell, multiplying the potential by 4irr 2 drp, the mass of the thin layer, and then 
integrating the resulting expression from r = to r = R, where R is the total radius of the spherical ball. This is 
given by, 



r R r2 _ 2 

U 2 = -2ttG1 2 sP / [Zn- -}4npr 2 dr 

Jo ( r + e ) - r 

After some algebraic manipulations, we obtain 



(35) 



„ ^,2 f R n R + r , R-r 
U 2 = -2nGl 2 p / [In + In 

Jo e 6 

, 2r + e, A 2 , 
—In \Airpr dr 



The total potential energy is given by, 

U = U 1 + U 2 



= -87r 2 Gl 2 p 2 [ [ r 2 dr(ln^-^ - ^LA) 



r dr(ln 



, R-r . 2r + e, 

In In , 

e e 

:In the e — > limit, it is not difficult to show that 



f l n A^ r ^ dr _ f i n ^L±l^ dr ^ (36) 

Je e Jo 6 

Therefore, the potential energy is given by, 

U= - 8n 2 Gl 2 p 2 [ I* r 2 dr(- 2(r ~ g) ) + f r 2 dr(ln^l + ln?—L)] (37) 
Je r Jo e e 



The last two integrals after change of variables amounts to evaluating standard table integrals of the form / x n lnx.We 
simply write down the results. 



" R R + r , 1 o R 2 „, . _ 5 



L 
I 



ln'——r z dr = -R 3 ln - + - R 3 ln2 R 3 

e 3 e 3 18 



R~ r 1 d3 , 11 n3 

In r 2 dr = -R 6 ln R 3 

e 3 e 18 
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Therefore, 



U = -8tt 2 G1 2 P 2 [ - R 3 In - + - R 3 In 2 - — R 3 }= -8tt 2 G1 2 sP 2 [ - R 3 In-- - R 3 In (5.16) ] 
3 e 3 9 J s L 3 e 3 

= -87r 2 GZ s y [ % R 3 ln{ U 



3 v 5.16e' 

With a redefined e, we can write it as, 



16 ^,2 2 n3 , R 3GM 2 l 2 R 
U = T^Gltp 2 R 3 In— = 7^r^-ln— 

3 s e R 3 e 



(38) 



If AT is the number of particles in the spherical ball and m is the mass of the elementary constituent particles, then 
M = Nm and we can write the total potential energy due to Randall-Sundrum brane-world induced correction term 
as 

rr rr 3GN 2 m 2 l 2 R 

U = U RS = 53— ^n- 39 

R 6 e 

The small parameter, e, in the formulae above should be interpreted as the minimum inter-particle distance in the 
spherically symmetric volume. 



VII. APPENDIX-II 



We want to calculate the extrema of the energy functionals considered in the main text of the paper. We will 
consider only the white dwarf stars. The results for neutron stars can be obtained simply by replacing the electron 
mass by the neutron mass. 

Most literature on white dwarf and neutron stars treat the ultra-relativistic limit of such stars taking the kinetic 
energy term as, K.E. = Nc\p\. We will first show that in the presence of brane corrections such systems do not have 
any stable quantum mechanical ground state. In this case the total energy including the brane world corrections as 
a function of R, is given by 

_, m AT , 3 Gm 2 n N 2 3Gm 2 n N 2 l 2 R 
E{R) = N c \p\ — - x 1 ^ '- x ln- 

Substituting p = hN J 3 , we obtain 



chN 4 / 3 3 Gm 2 n N 2 3Gm 2 n N 2 l 2 , R 

= — — 5 X — i W^ xln 7 



With N max and a, defined as 



5 ch 

N max - X ^2 ! a = ( N / N ma X ) 2/3 

we can write the total energy given above as 

,^/^\ ,..,//?,! a 15aZ, , R, 
E(R) = chN^ 3 [- --- x In-} 

The extremum is obtained by equating the the derivative of E(R) with respect to R to zero with N held fixed. 

dE chN A / 3 45o^ , R 

dR^-^^ {l - a -^ xln 7[ ) = {] 

where t\ — je ; 7 being some positive number, 1 < 7 < 10 . To find the nature of the extremum, we calculate the 
second derivative of E(R) with respect to R. 

d 2 E 2chN A / 3 , _ 45ai? , R, cTiN 4 / 3 ,90 a I 2 , R, 

m = —R^ {1 - a -^ xln 7?--R^ { -® ±xln ^ 
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where 62 = 7161 ; 71 being some positive number, 1 < 71 < 10 . The quantities e, e\ and €2 as an order of 
magnitude can be taken to be equal to j^js- In the equation above, the first term is zero from the condition that 
^| = 0. The second term is strictly negative. Therefore, the extremum of the function, E(R), is a maximum. The 
energy E(R) as a function of R, starts from —00 at R = 0, rises very sharply to a positive value and develops a 
maximum very close to R = 0, and then from the positive side of energy it slowly approaches zero at R = 00. 
Next we consider the total energy functional given by 

3 Cm 2 /V 2 

E(R) = K N + U N (R) = Nc(p 2 + mlc?) 1 ' 2 - Nm e c 2 - ± x - " 

IX 



Substituting p - 



hN 



1/3 



R ' 



we obtain 



E(R) = N c ( R2 + rale 2 ) 1 ' 2 - Nm e c 2 - - x 1 



The extremum with respect to R is given by the condition, 



dE 
dR 



zUN 4 ' 3 



R 2 



N l ' 3 h 



(h 2 N 2 / 3 + mlc 2 R 2 y/ 2 5 



3 GmlN 2 ^ 
x 



cfi 



In terms of the parameters, Rq ( 



hN 1 ' 3 
m e c 

dE_ 
dR 



) and a , the equation above simplifies to the form 



R 2 



(l + f) 1/2 



To find the nature of the extremum, we calculate the second derivative of E(R) with respect to R. It is given by, 



d 2 E _ 2chN 4 ^ 
dR 2 ~ ~ 



(l + f) 1/2 



~~R 2 ~ 



-R 



i?o 2 (l + f) 3/2 



The first term is zero because of the condition, |^ = 0. Therefore, 



d 2 E 
dR 2 



2hN 4 / 3 



R 



i?o 2 (l + f) 3/2 



which is strictly positive implying that the extremum is a minimum. The energy E(R) as a function of R, starts from 
oo at R — and rolls down to a minimum at a negative value of E(R) and then from the negative side of energy it 
slowly approaches zero at R = oo. 

Now we consider extremum of the total energy functional given by, 



E(R) = K N + U N (R) + U RS = Nc (p 2 + mlc 2 ) 1 ' 2 - Nm e c 2 - | x _ 



•ZGls 2 mlN 2 , R 
- x In— 



Substituting p - 



hN 



1/3 



we obtain 



h 2 N 2/3 

E(R) = N c ( R2 + m 2 e c 2 ) 1/2 - Nm e c 2 



3 GmlN 2 

- x 

5 R 



3Gls 2 m 2 n N 2 
R 3 



; R 

x In— 

e 



The necessary condition for the extremum is given by. 
dE(R) _ c N b ' 3 h 2 



dR 



i?3(Z H^ +m 2 c 2)l/2 5 



3 GmlN 2 9Gls 2 m 2 n N 2 , R 3Gls 2 mlN 2 
x ^ 1 —^ x In— 



R 2 



R 4 



R 4 



The last term, in the equation above, can be absorbed in the third term containing the logarithm by simply redefining 
the parameter e. As before, we denote it by e\. With the redefined e, the equation above takes the following form, 



dE_ 
dR 



zUN 4 ' 3 
R 2 



N^h 3 Gm 2 n N 2 / 3 §Gls 2 m 2 n N 2 l 3 , R 

— x - — x In 

(h 2 N 2 / 3 +m 2 c 2 R 2 y/ 2 5 c% R 2 ch e 1 
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In terms of the parameters, Rq and a , the equation above simplifies to the form 



dE 
dR 



R? 



1 



15a I 2 , R 

55 — ; a t; — x In — 

(l + f)V 2 R 2 



In the main text of the paper, it is shown that 

solves the equation above. The parameter (3 is defined in the text as, 

„ 15/2 R \U 2 S , Arl/3 
(3 = — f x In— w — f x /^iV 1/3 

For l s = 1mm and N « 10 57 , we find that (3 ~ Therefore, the solution of the equation for the cxtrcmum ean also 
be written as 



To make the brane world parameter dependence explicit, we can also write it as 

R 2 a 2 { R 2 > 



(40) 



(41) 



To find the nature of the extremum, we need to calculate the second derivative with respect to R of the function 
E(R). This is given by, 



d 2 E _ 2chN 4 / 3 

ZR2 - ^3 



1 



15a l 2 s , R 

59 — ; a n — x in — 

(1 + ^)V2 R 2 e 1 J 



chN^ 
R 2 



-R 



30a I 2 R 



The first term in the equation above is zero because of the condition, = 0. With l s = 1mm = 0.1cm and e 2 = jjr/j, 
we can write the equation above as, 



d 2 E 
dR 2 



chN 4 / 3 
R~ 



0.3a InN 1 / 3 



^(l + f) 3/2 



R 4 



Using Eq.( ), this can be written as, 

d 2 E chN 4 / 3 a 3 



dR 2 



R 



0.3a 2 InN 1 / 3 



Rt + 0{P s ) Ri(a 2 -l) 2 + O(l 2 R 2 ) 



(42) 



For all admissible values of the parameters for white dwarfs as well as neutron stars, the expression above is strictly 
positive. Therefore, the extremum is a minimum. Note that the equation for extrema is a cubic equation in R 2 , and 
therefore, there should be two more solutions. These two solutions up to order j3 are real and coincident. These two 
solutions are given by, 



R 2 



= 0(P) = 0{^) 

-fin tic. 



(43) 



These solutions correspond to a maximum of the total energy E(R). The energy E{R) as a function of R, starts from 
— oo at R — 0, rises very sharply to a positive value and develops a maximum very close to R — 0, it rolls down to a 
minimum at a negative value of E(R) and then from the negative side of energy it slowly approaches zero at R = oo. 



[1] L.Randall and R.Sundrum, Phys.Rev.Lett. 83, 3370 (1999) 



15 



[2] L.Randall and R.Sundrum, Phys.Rev.Lett. 83,4690 (1999) 

[3] T. Shiromizu, K. Maeda and M. Sasaki, Phys. Rev. D62,024012 (2000); M. Sasaki, T. Shiromizu and K. Maeda, Phys. 
Rev. D62,024008 (2000); Roy Maartens, Phys. Rev. D 62, (2000) 084023; J.M. Cline, C. Grojean and G. Servant, Phys. 
Rev. Lett. 83 4245 (1999); N. Kaloper, Phys. Rev. D 60 123506 (1999); P. Binetruy, C. Deffayet, U. Ellwanger, and D. 
Langlois, Phys. Lett. B477, 285 (2000); E. E. Flanagan, S.-H. Tye, and I. Wasserman, Phys. Rev. D 62 044039 (2000); 
S. W. Hawking, T. Hertog, and H. Reall, Phys. Rev. D 62 043501 (2000); S. Nojiri, S. D. Odintsov, and S. Zerbini, Phys. 
Rev. D 62 064006 (2000); S. Nojiri and S. Odintsov, Phys. Lett. B 484, 119 (2000); L. Anchordoqui, C. Nunez, and K. 
Olsen, J. High Energy Phys. Phys. 10, 050 (2000); S. Nojiri and S. Odintsov, Phys. Lett. B 494, 135 (2000); S. Gubser, 
Phys. Rev. D 63, 084017 (2001); T. Shiromizu and D. Ida, Phys. Rev. D 64, 044015 (2001). 

[4] Roy Maartens, D. Wands, B.A. Bassett, and I.P.C. Heard, Phys. Rev. D 62, 041301 (2000); E.J. Copeland, A.R. 
Liddle, and J.E. Lidsey, |astr o-ph/0006421 E.J. Copeland, A. R. Liddle and J. E. Lidsey, Phys. Rev. D 64, 023509 
(2001) |astro-ph /0006421|; G. Huey and Lidsey, Phys. Lett. B 514,217(2001); V. Sahni, M. Sami and T. Souradeep, 
Phys. Rev. D65,023518(2002); A. S. Majumdar, Phys.Rev. D64 (2001) 083503; J. W. F. Valle, As tropart. Phys. 18, 
287(2002) astro -ph/0111417| ; K. Dimopoulos, Nucl. Phys. Proc. Suppl.95, 70(2001) astro-ph/0012298 |; M. Sami, N. D ad- 
hich and Tetsuya Shiromizu, Phys. Lett B 68, 118(2003) hep-th/0304187 ; M. Sami and V. Sahni, hep-th/0402086 M. 
Sami and N. Dadhich, hep-th/0405016 Shinji Tsujikawa and A. R. Liddle, JCAP 0403 (2004) 001 astro-ph/0312162 . 

[5] C.Germani and R.Maartens, Phys. Rev. D64 124010(2001) 

[6] R.Maartens, Brane-World Gravity, gr-qc/0312059 

[7] T.Wiseman, Phys. Rev. D65 124007(200 2); |hep-th/0111057| 

[8] T.Damour, Gravitation and Experiment, gr-qc/9711061 

[9] E.H.Lieb et. al. The stability of matter:From atom to stars, Selecta of E.H.Lieb, Ed. W.E.Thirring, Springer- Verlag, 1991 
[10] S.L.Shapiro and S.A.Teukolsky, Blacl Holes, White Dwarfs, and Neutron Stars; John Wiley and Sons, Inc, 1983 
[11] J.Hartle, Physics Reports 46C 200(1978) 



